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Abstract We describe a novel high aspect ratio ra- 
diofrequency linear ion trap geometry that is amenable 
to modern microfabrication techniques. The ion trap elec- 
trode structure consists of a pair of stacked conducting 
cantilevers resulting in confining fields that take the form 
of fringe fields from parallel plate capacitors. The confin- 
ing potentials are modeled both analytically and numer- 
ically. This ion trap geometry may form the basis for 
large scale quantum computers or parallel quadrupole 
mass spectrometers. 

PACS: 39.25.-fk, 03.67.Lx, 07.75. +h, 07.10+Cm 



1 Introduction 

The ion trap has become an essential tool in several 
areas of physical science, including mass spectroscopy 
, atomic frequency standards [2] , precision atomic and 
molecular measurements '3' , studies of fundamental quan- 
tum dynamics W and quantum information science 0E1 ■ 
Many of these applications would benefit from miniatur- 
ized and multiplexed ion trap electrode structures well 
below the typical millimeter to centimeter scale. Further- 
more, smaller electrode dimensions offer the potential for 
stronger confining forces. 

In this paper, we consider theoretically the electri- 
cal characteristics of a new type of micrometer-scale ra- 
diofrequency (RF) Paul ion trap fabricated using semi- 
conductor micromaching and lithographic techniques such 
as micro-electro-mechanical-systems (MEMS) and molec- 
ular beam epitaxy (MBE). Such a device may enable new 
applications of ion trap technology such as "quantum 
CCD" scalable quantum computers optical cavity- 
QED with a localized single atom and mul- 

tiplexed quadrupole mass spectrometers that could be 
orders of magnitude smaller than previous devices |12| . 

There has been much recent progress in the miniatur- 
ization of neutral atom electromagnetic trapping struc- 
tures, involving, for example, micrometer-scale current- 



carrying wires on a substrate resulting in Bose-Einstein 
condensates on a microchip |13j . Microscopic ion trap 
electrodes present their own challenges, as the confining 
forces are orders of magnitude stronger than those for 
neutral atom traps. Consequently, such ion traps will re- 
quire greater control of unwanted or noisy electrode po- 
tentials, including the presence of thermal electric fields 
|14lll5j . residual charge on exposed insulating barriers, 
and "patch" potentials from inhomogeneities on the elec- 
trode surfaces |UIITB] . None of these potential pitfalls ap- 
pears fundamental, and such problems will only be over- 
come by testing various materials and approaches. We 
focus here on novel features of a proposed high aspect- 
ratio ion trap geometry and the resulting confining po- 
tentials. 



The physical parameters of a model of the linear mi- 
crotrap are discussed in Sec. |21 along with a discussion of 
design considerations and issues with heating and power 
dissipation in semiconductor materials. Section |31 con- 
tains a discussion of the RF ponderomotive potential of 
the linear microtrap model with results from numerical 
simulations of the potential. A geometrical efficiency fac- 
tor is calculated, showing the performance of the linear 
microtrap as compared to an ideal quadrupole potential. 
The static potential used for axial confinement in a linear 
trap is discussed in Sec.^along with results from numer- 
ical simulations and comparison to an ideal hyperbolic 
trap. The total potential along with examples of how to 
use the various geometric efficiency factors to calculate 
the trap frequencies of a given geometry are given in Sec. 
121 The principal axes of the linear microtrap, which de- 
termine the efficiency of laser cooling ions in the linear 
microtrap, are evaluated in Sec. El A method for rotating 
the axes for more efficient cooling is given. 
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2 Model Description 

2.1 Basic Model 

The design of this new type of micrometer-scale RF trap 
is constrained by conventional semiconductor fabrica- 
tion techniques, the need for clear laser optical access, 
and the characteristics of electrodes needed for linear 
traps. The design, illustrated in Fig. ^ is a two-layer 
planar geometry where both layers are divided into sep- 
arate electrodes. The division of each layer into six elec- 
trodes accommodates both the RF potentials and the 
static potentials needed to create a linear Paul trap 17j . 
This planar design is compatible with conventional pho- 
tolithography techniques to define the electrode pattern. 
Each electrode is a cantilever anchored to an electrically 
isolated, conductive substrate and suspended from both 
sides of the planar structure. This ensures that there are 
no insulators near the center of the trap that could ac- 
cumulate uncontrolled charge. Ions will be trapped in 
the space between the tips of each cantilever, along the 
2-axis in Fig.^ near the center of the middle electrode. 



the model is h. In order to electrically insulate the center 
from the end-cap cantilevers, a small gap is introduced of 
width g. This allows for separate potentials to be apphed 
to all twelve cantilevers, or electrodes. Static voltages 
are applied to both layers on the four end-cap electrodes 
on either side of the center cantilevers to provide axial 
confinement, as shown in Fig.|2Ib). 

The potentials of the LMT can be separated into 
two parts for analysis. The first part is the ponderomo- 
tive potential generated by the RF voltages. In the limit 
where gap width g is much smaller than a, b and c (Fig 
|2Ib)), the RF potential is approximately independent of 
z near the center of the trap. In the cross-sectional plane 
at 2 = 0, this RF potential generates a two-dimensional 
trapping pseudopotential and is discussed in Section |21 
The second part is the potential generated by applying 
static voltages to the end-cap electrodes. This potential 
provides axial confinement for ions in the center of the 
trap and is described in Section^] Note that the end-cap 
electrodes have both the RF voltages applied to reduce 
the z dependence of the RF field near the center of the 
trap and static voltages to create the end-caps. The cen- 
ter electrodes are all assumed to be held at static ground. 




Fig. 1 A three-dimensional drawing of the linear microtrap. 
A string of ions would lie along the z-axis as drawn. 



The cross-section of this linear microtrap (LMT) model 
at the center of the trap {z = 0) is shown in Fig. [2Ja) . 
The thickness of each layer is labeled w; the layer sepa- 
ration is d] the tip-to-tip separation of the cantilevers is 
a. Two ratios are useful for characterizing the behavior 
of the electric potentials: the trap aspect ratio, or the 
ratio of the tip-to-tip cantilever separation to the layer 
separation a = a/d, and the ratio of the layer separation 
to the layer thickness 6 = d/w. An RF voltage is applied 
between each set of diagonally opposing electrodes as 
shown in Fig. [S^a). 

A top view of the linear microtrap model is shown 
in Fig. EJb) . The width of the center cantilevers along 
the 2;-axis of the trap is labeled b; the width of the end- 
cap cantilevers is c; and the length of the cantilevers in 
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Fig. 2 (a) A schematic diagram of the linear microtrap de- 
sign showing the side view. The dimensions are labeled as are 
the RF voltages applied to the electrodes, (b) The top view 
of the linear microtrap with dimensions and static voltages 
as shown. 



2.2 Fabrication Considerations 

The linear microtrap model is designed to simulate a 
trap design that can be fabricated using conventional 
micro-processing techniques. The sizes of the electrode 
features that will be analyzed in this model are typical of 
current fabrication processes. There are several different 
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processes that could be used to fabricate these micro- 
traps: a) Silicon-based microelectromechanical machin- 
ing (MEMS) techniques; b) GaUium-Arsenide (or other 
suitable material) based molecular-beam epitaxy (MBE) 
grown wafers and associated etching processes; or c) 
other relevant techniques such as anodic wafer bonding 
or flip-chip technologies. The length of the cantilevers is 
limited by allowable mechanical vibrations in the can- 
tilevers themselves, as well as limits to the mechani- 
cal stabihty of the cantilevers under electromechanical 
forces due to the applied RF and static voltages. The 
mechanical forces exerted on the cantilevers can be ap- 
proximated using structural cantilever analysis ^H] . Fol- 
lowing this analysis, the spring constant of the center 
rectangular cantilever can, for example, be expressed as 



(1) 



where E is the Young's Modulus of the relevant material. 
The force on one cantilever due to an applied potential 
difference Vq between layers can be approximated as the 
gradient of the potential in a parallel plate capacitor of 
area A = hb and plate separation d. 

^U(Japacitor 

" dd 

^ £0 a_ fhwi 

2 dd\ d 
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(2) 



Although the actual force is distributed across the length 
of the capacitor, by approximating the force as being 
concentrated at the tip, one can find an upper bound on 
the cantilever tip deflection. Treating the cantilever as 
a classical spring with the force applied at the tip, and 
using the spring constant from Eq. ^ the maximum tip 
deflection x^?^ can be approximated as 



2£ofeVo' 
Ed'^w^ 



(3) 



A typical deflection for a GaAs cantilever with E = 
85.5GPa and dimensions h — 100/im, d — 2/im, w = 
2/im, with an applied voltage difference of Vq = 20V, is 
x^j^^ = 260nm. The resonant frequency of the cantilever 
can also be calculated 18 as a function of the material 
density p, the Young's Modulus, the cantilever width w 
and the length h: 



w„6/27r = 0.162/E7pT7 



(4) 



which for GaAs {p — S.Slgm/cm'^) is ujyit,/2TT « 130kHz 
for the same dimensions as previously discussed. 

For an RF potential Vq cos{f^Tt) applied to the can- 
tilever electrodes, the amplitude of the tip deflection in 
Eq. Olis expected to be further reduced by a Lorentzian 
factor of ojl^jj/n^ ^ 1. Here, it is assumed that the RF 



frequency is far from resonance, or Qt ^ ^vib/Q, where 
Q ^ 1 is the quality factor of the mechanical resonance 
[T^. While the above electromechanical forces do not ap- 
pear troublesome, the actual forces may be considerably 
higher due to free charges on the electrode layers that 
are driven by the applied potentials. In any case, it may 
be necessary to isolate the cantilevered electrodes from 
noisy electrical signals near the mechanical resonance. 

The trap strength may be limited by the maximum 
voltage that can be applied to the electrodes before the 
occurrence of electric field break-down. The theoreti- 
cal limit to the breakdown voltage is dependent on the 
bandgap of the semiconductor material and, for Si and 
GaAs, is on the order of 40-50 V/pm |5n| and for silicon 
nitride, on the order of 300 V/ pm For a layer sepa- 
ration of 2pm, the maximum applied voltage is expected 
to be of order Vq = lOOV. 



2.3 RF Dissipation and Thermal Fields 

The fabrication considerations for the implementation 
of this new type of linear microtrap suggest that highly 
doped semiconductors could be used as electrodes. Be- 
cause doped semiconductors have a resistivity several 
orders of magnitude greater than the metal conductors 
typically used in ion traps, it is necessary to estimate 
the power dissipation of the microtrap due to RF losses 
in the cantilevers. Additionally, the finite conductivity 
of semiconductor materials will lead to thermal electric 
fields that will generate heating of the quantized motion 
of ions in the center of the trap. 

The RF dissipation can be estimated with a sim- 
ple model of lumped circuit elements, since the trap 
structure is much smaller then the RF wavelength. Each 
RF electrode is modeled as a small series resistance R 
shunted by a capacitance C at the trap; inductance of the 
electrodes is assumed negligible compared to 1/(CJ7|.). 
In addition, RF loss in the insulator separating the elec- 
trodes contributes to a parallel resistance characterized 
by the loss tangent tan (5. Assuming RCnT,ta.nS <C 1, 
the power loss is 



V^Cf2T 



{RCf2T + tanS). 



(5) 



For values envisioned here, Vq ^ 20V at J7r/27r ^ 50MHz, 
C - lOpF, tan 5 ~ 0.0002 and R ~ 1012, resulting in a 
power dissipation of Pd ^ 40mW per electrode. 

Additionally, Johnson noise in the electrodes will gen- 
erate thermal electric fields that will cause heating of 
the quantized ion motion. A simple model can be used 
to calculate the heating due to the resistivity of the trap 
electrodes |Hlll5| . For an ion held at a distance z from a 
conductive plane, the heating rate is given by 

dE 
'dt 



hujn 
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where u is the secular frequency and n is the average vi- 
brational quantum number of an ion in the trap. In the 
limit where the conductor thickness w is much smaller 
than the distance to the ion z, and both dimensions are 
smaller then the skin depth 5 of the conductor {w <^ 
z <^ 5), the resistance R in Eq. Elis frequency indepen- 
dent: R ~ pz/{zw), where p is the material resistivity. 
Here, the effective volume of the conductor contributing 
to the thermal fields is of order z^w. Again, using typical 
values for doped semiconductors, the skin depth J is a 
few hundred micrometers, the thickness of the conduc- 
tor is 2/im and the ion is 20/im from the conductor. In 
this limit, using a secular frequency of u;s/27r = lOMHz, 
and ^^^Cd+ ions, Eq. El predicts a thermal heating rate 
of about 10 quanta/sec. Since this model pertains to 
fluctuating uniform thermal electric fields from a single 
conducting plane, the actual thermal electric fields are 
expected to be much smaller because the trap structure 
surrounds the ion with a high degree of symmetry, re- 
sulting in some degree of cancellation of thermal fields 
from opposite electrodes. In any case, the heating rate 
will likely be limited in practice by fluctuating patch 
fields on the electrode surfaces [T^. 

3 RF Ponderomotive Potentials 



the RF potential is approximately uniform along the z- 
axis near the center of the trap, it can be described in 
the z — Q plane as a function only of x and y, reducing 
the calculation of the RF potential to two dimensions. 



3.2 Hyperbolic Electrode Model 




Fig. 3 The hyperbolic electrode geometry is used as a basis 
for comparing the linear microtrap. The characteristic dimen- 
sion of the hyperbolic electrode geometry is the radius _Ro as 
shown. 



3.1 Time- dependent RF potentials 

As described above, the analysis of the potentials in a 
linear RF Paul trap can be divided into the transverse 
RF trap generated by RF voltages applied to the appro- 
priate electrodes, and the axial trap and transverse anti- 
trap generated by static voltages applied to the end-cap 
electrodes. Focusing first on the time-varying potential 
generated by the RF voltages, the analysis can be sim- 
plified by using a pseudopotential approximation. The 
motion of an ion in an RF potential of the form 



<?(x, y, z, t) = V{x, y, z) cos{f2Tt) 



(7) 



can be approximated using a ponderomotive pseudopo- 
tential [SI: 



'WV{x,y,z)\^ 



(8) 



Ion motion is in the pseudopotential can be approxi- 
mated as secular harmonic motion |S] with frequency 



\^V{x,y,z)i 



(9) 



The micromotion due to the time dependence of the RF 
potential is small in the limit where q = 2\/2LOp/ Qt <C 1 

m 

Since the secular ion motion is dependent only on 
the gradient of V{x,y, z), it is possible to calculate the 
effective (or ponderomotive) potential of the linear mi- 
crotrap using an electrostatic analysis. Moreover, since 



One common configuration of a linear Paul trap con- 
sists of four infinitely long hyperbolic electrodes. This 
hyperbolic electrode model will be used as a standard of 
comparison for the linear microtrap. The cross-section 
of hyperbolic electrodes with a characteristic radius Rq 
is shown in Fig. |31 For the potentials applied according 
to Fig. 13 the exact potential amplitude is 



Vi,yp{x',y') 



2Rl 

Vor^ 
2 i?g 



(10) 



cos 261'. 



where the coordinate system {x',y') is indicated in Fig. 

m 

The pseudopotential that corresponds to this hyper- 
bolic potential is calculated using Eq. |S| 



^hyp 



° {x" + y'') (11) 



The secular frequency of a ion moving in this pondero- 
motive pseudopotential is therefore 



^p,hyp 



eVo 



V2m[2tRl 
/ eVoq 



AmRl' 



(12) 
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3.3 Linear Microtrap Transverse Potential Analysis 

The microtrap potential amplitude V^mt computed 
near the center of the trap. This potential is then de- 
composed as an infinite set of cylindrical harmonics |^ : 



V,„,{r,e') = Vu 



^ Cm {r/ro)"" cos{m0') 



+ J2Sn{r/rorsm{n0') 



(13) 



where Cm and Sn are expansion coefhcients and Q' is 
taken as the angle from the x' axis. The characteristic 
radius over which the potential is approximated by this 
expansion is tq. 

The C2 coefficient provides a comparison between the 
potential of the linear microtrap and the quadrupole po- 
tential of the hyperbolic electrode geometry of radius 
tq. Other nonzero coefficients in the expansion of Eq. 
1131 describe the anharmonic character of the microtrap 
potential. Symmetry considerations reduce the number 
of terms allowed in the expansion. Given the potential 
amplitude of ±Vb/2 applied to opposite electrodes as 
shown in Fig.[21[a), the potential is antisymmetric along 
the lines x = and y = and symmetric in reflection 
about the origin leading to the only non-zero terms in 
Eq.HSas m = 2,6,10, ... and n = 4, 8, 12, . . .. 

The expansion coefficients are calculated by numeri- 
cally evaluating the LMT potential using finite element 
analysis or other appropriate numerical field simulators 
and calculating the overlap integrals within a circle of 
radius of the potential V^mt with the cylindrical har- 
monics (r/ro)™ cos(m6'') and (r/ro)" sin(n0') |23j . 

A geometric efhciency factor r\ can be used to com- 
pare the microtrap potential with the quadrupole poten- 
tial of the hyperbolic electrodes of comparable size. The 
size of the linear microtrap is given by the distance from 
the center of the trap to the nearest point on the tip of 
the electrodes 4ff = 1/(0/2)^ ^ {d/2Y. Then, 77 is de- 
fined as the ratio of the quadrupole part of the potential 
generated by the LMT V^'^^ and a hyperbolic trap with 



2C2e 



The quadrupole portion of the linear microtrap can there- 
fore be written in a form differing from the hyperbolic 
electrode potential fEa. 111(1 by only the geometric factor 

V'^) (15) 



LMT ' » / " ""^ ' 



2P 



The ponderomotive potential for the microtrap can then 
be evaluated using Eq. |S1 



2^2 



e V^ri 



y 



(16) 



Finally, the effective secular frequency of an ion in the 
linear microtrap is only modified by the factor 77 from 
the form of the secular frequency in the trap due to the 
hyperbolic electrodes (Eq. 112(1 . With this form of the 
secular frequency, one can compare the trap strength 
and performance of the linear microtrap. 



eVori 



P,LMT 



V2mf2Til 



(17) 



off 



The equipotential lines of the calculated ponderomo- 
tive potential are shown in Fig. ^ along with the po- 
tential magnitude indicated by a gray-scale. Note that, 
although the cantilever geometry does not have cylin- 
drical symmetry, the pseudopotential is approximately 
circular within a distance on the order of one-eighth the 
tip-to-tip separation a as will be shown from the nu- 
merical results in Sec. 13.41 where C2 is found to be the 
dominant term in the expansion at this distance from 
the center. 




Fig. 4 Equipotential lines of the pseudopotential Vlmt 
the z — plane for aspect ratio a = 10 and ratio of layer 
separation to layer thickness of 5 — 1. The ponderomotive 
potential reaches a maximum along the y axis near l^s- The 
contour lines are spaced on a linear scale and are shown to 
illustrate the circular nature of the ponderomotive potential 
at the center of the trap. The gray-scale shading is also on a 
linear scale. 



3.4 Finite Element Analysis Method 

The class of finite element analysis solvers that is used 
here divides a two-dimensional space into a series of tri- 
angles to calculate the linear microtrap potential. The 
two-dimensional finite element analysis package in Mat- 
lab version 6.5 was used to calculate the RF potentials. 
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The results were compared with the two-dimensional 
projection of potentials calculated using two different 
three-dimensional finite element analysis packages, Maxwell 
3D from Ansoft, and Opera 3D from VectorFields, and 
found consistent. The field is approximated at each ver- 
tex on the triangles, then an interpolation is made within 
each element to calculate the field on an rectangular 
grid. Different trap configurations are analyzed using the 
method described above and the ratio rj of the micro- 
trap potential to the quadrupole hyperbolic potential is 
shown in Fig. evaluated at a radius of rp = a/8. The 
uncertainty of the simulation data is less then 5% and 
is due primarily to a finite grid spacing and the finite 
bounding box size. The solid line in the figure is an an- 
alytic solution for cantilevers of infinitesimal thickness. 

One can see that as the trap aspect ratio a = a/d 
increases, the geometric factor rj approaches a constant, 
non-zero value. The asymptotic value can be evaluated 
using complex analysis techniques and is described in 
Appendix ^ The result from Eq. 1511 for large a is rj — 
I/tt. Additionally, as the aspect ratio approaches one, 
the trap becomes more like the hyperbolic electrode ge- 
ometry. The other degree of freedom of the linear mi- 
crotrap is the ratio of the layer separation to the layer 
thickness, 6 = d/w. Note that the strength of the mi- 
crotrap decreases as the layer thickness decreases with 
respect to the layer separation. 




40 60 
Aspect Ratio a 



Fig. 5 The ponderomotive potential geometric efficiency 
factor ry as a function of the ratio of the tip-to-tip separa- 
tion to the layer separation: the aspect ratio a. The other 
degree of freedom is the ratio of the electrode separation to 
the layer thickness, 5 = d/w. The solid line is an analytic 
solution for rj found using complex analysis techniques with 
5 — > oo and is valid for a S> 1 



The higher-order coefficients of the expansion shown 
in Eq. E| for the potential Vlmt ^^'^ shown in Fig. El 
The dominant higher-order term is 54, which, at a fixed 
radius of tq = a/8, is only a few percent of C2. The two 
next largest terms are also shown although the magni- 
tude is small enough to be negligible when considering 
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Fig. 6 The two largest higher-order terms of the expansion 
in Eq. 1131 shown as a ratio over C2 for various trap aspect 
ratios a — a/d and given as a function of the layer separation 
over the layer thickness 5 — d/w evaluated at ro = a/8 



ion motion. The relationship between the C2 and the 
next three largest terms of the expansion as a function 
of the aspect ratio a and 6 is shown in Fig |H| Coefh- 
cients S4, Cg appear to approach an asymptotic value as 
the trap aspect ratio increases. The ratios of all higher- 
order terms to the coefficient C2 (Cm/C2 and S'„/C2) 
for TO,n > 6 are less than 10"'^. 

The absolute depth of the ponderomotive RF trap 
is also of interest when considering ion loading and col- 
lisions with background gas. The trap depth is defined 
as the maximum height of the ponderomotive potential 
barrier along the weak axis of the trap and is plotted 
in Fig. [7| A trap frequency of I^t/Ztt = 50MHz and 
the mass of ^^^Cd"*" were used to calculate the depth, 
given in scaled units of [K-/im^/V^]. To find the depth 
of a specific trap, the data must be multiplied by the 
applied voltage Vq in [V^] and divided by the square of 
the absolute tip-to-tip separation in [/.im^]. The depth 
asymptotically approaches a constant value of approxi- 
mately 2400K-/im^/V^ for large cross-sectional aspect 
ratio as can be found from the analytic solution (Eq. 
154(1 . The size of the ponderomotive trap rmax is char- 
acterized by either the distance of the maximum in the 
ponderomotive potential from the center of the trap or 
a/2, whichever is smaller. As the trap aspect ratio in- 
creases Tmax is determined by the maximum in the RF 
pseudopotential along the y-axis and is approximately 
half the tip-to-tip electrode separation 0.5a. 

Since the ponderomotive potential within the region 
r < rmax will trap ions, the expansion of the potential 
from Eg. 1131 within that entire area is also of interest. 
The expansion of the potential within a circle of radius 
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Fig. 7 The scaled trap depth as a fiinction of the trap as- 
pect ratio a = a/d and the ratio of layer separation to layer 
thickness 5 = d/w. The trap depth is scaled to the tip-to-tip 
separation a, in micrometers and to the applied voltage Vb- 
The analytic result is shown as a solid line with 5 —> oo and 
is valid for a S> 1. 



'''0 = ''max contains a larger contribution from the higher- 
order coefficients than an expansion fixed at ro = a/8 as 
iUustrated in Fig.|Sl The higher-order coefficients for the 
expansion of the hnear microtrap potential are shown in 
Fig. evaluated at ro = ^max- 
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Fig. 8 The dependence on the expansion coefficients Cm 
and Sn as a function of ro/(a/2). The higher order terms 
become significant as the overlap integrals cover more of the 
area between the electrodes. The geometry used was a — 20 
and (5 = 1, a worst case scenario from Fig. |5] The dashed 
vertical bar indicates ro = a/8. 



3.5 Residual Axial Ponderomotive Potential 

The previous analysis is based on the assumption that 
the linear microtrap electrodes are infinitely long in the 



Fig. 9 The three largest higher-order terms of the expansion 
in Eg. 1131 evaluated within a radius ro — rmax, where rmax 
is the maximum of the ponderomotive potential. The coeffi- 
cients are shown as a percentage of the largest term C'2 for 
various trap aspect ratios a = a/d and given as a function 
of the layer separation over the layer thickness S — d/w with 
an error of 5%. 

z-dimension. However, the actual trap has finite elec- 
trode lengths, labeled b and c in Fig. |2Ib), which to- 
gether with the small electrode gaps (labeled g in Fig. 

lead to a small ponderomotive potential in the z direc- 
tion. The magnitude of this axial ponderomotive poten- 
tial can be compared to the transverse ponderomotive 
potential V'lmt of Eq. 1161 To find the axial contribution, 
the entire three-dimensional RF potential V^^,j,{x,y, z) 
must be computed. Once found, one can use the pon- 
deromotive potential approximation Eq. [S] to calculate 
the trap frequency along the z-axis. 

The gradient of the three-dimensional potential is 
found, then the pseudopotential is evaluated. A Taylor 
expansion of the pseudopotential along the z axis (about 
z = 0) gives the coefficient for the harmonic term in 
the ponderomotive potential: 

H.^^-^{\VV,^,{x,y,zfy (18) 

The details of the three-dimensional potential calcula- 
tion are given below, but the method is similar to the 
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two-dimensional finite difference analysis. Typically the 
data is extracted along the z axis and then fit to a 
quadratic polynomial to find the coefficient Hz ■ This co- 
efficient allows one to make a comparison between the 
quadrupole trapping pseudopotential in the z — cross- 
sectional plane, and the ponderomotive potential along 
the z-eods. This three-dimensional ponderomotive poten- 
tial is similar to the transverse potential of Eq. 1161 with 
the addition of the term. 



Amt ix,y,z) 



1 err 



(19) 



where CTz = Hzt-'^^/rf' is the ratio of the residual ax- 
ial to transverse ponderomotive potential. The resulting 
frequency along the z-axis is uJz = \/^^p_-LmT 

The results from the numerical simulation in Fig. 1101 
are given for a cross-sectional aspect ratio of a = 20 
and for 5 — 1 (ratio of the layer separation to the layer 
width) . The ponderomotive potential along the z-axis is 
shown in Fig. ^2 to illustrate the degree to which the 
notch gap g contributes to the residual potential at the 
center of the trap. Since ct^ ^ 1, the ponderomotive 
contribution to the potential along the z-axis can be ne- 
glected. 



-1e-15 



-1e-14 



-1e-13 



-1e-12 

1 2 3 4 5 5 7 8 
Center electrode length b/f^^ 

Fig. 10 The ratio of the residual axial frequency to the 
transverse ponderomotive frequency as a function of the 
center electrode length. The end-cap electrodes were fixed at 
5£eff with a fixed gap spacing of 1/lO^off, ot = 20, and 5 = 1. 



4 Static Potentials 

4--1 Hyperbolic Geometry 

Like the two-dimensional potential in Sec. 13.21 the static 
potential used to confine the ions along the z-axis in the 
linear microtrap can be compared to a three-dimensional 
idealized hyperbolic electrode potential. FigurelT^shows 
an elliptical hyperbolic electrode geometry where a;o , yo, 
and zq are the distances along the principal axes of the 
ellipse from the center of the trap to the electrodes. The 




T 1 1 1 r 


/ 


b=1S0nm 








b=lOOum \ 










: 1. 


b=40nm \ 



z [nm] 

Fig. 11 Illustration of the change in the residual axial pon- 
deromotive potential for various center electrode lengths (6). 
The potential along the z-axis is shown for various center 
electrode lengths where the end-cap electrodes have been 
fixed at 100/im. 




Fig. 12 Three-dimensional hyperbolic electrodes are shown 
here. The electrodes along the z-axis are held at a voltage 
of C/o, while the center electrode is grounded. The potential 
has an elliptical cross-section in the xy plane corresponding 
to e = 0.86 and, for C/o > 0, is trapping along the z-axis, but 
anti-trapping along x and y, valid for < e < 1. 



potential within the electrodes, up to a constant term, 
is 



t^hyp = {-ex^ - (1 - e)y^ 



(20) 



where — Zq+ exQ and ex^ = (1 — e)2/o- The geometric 
anisotropy factor e is related to the eccentricity of various 
conic sections that can be superimposed on the three- 
dimensional hyperbolic electrode structure. The special 
case where e = 1/2 corresponds to circular symmetry 
in the xy plane. For values of < e < 1 and Uq > 0, 
the potential is trapping in z and anti-trapping in the 
xy plane, as shown in the figure for e = 0.86. Outside 
of that range, the axes in the figure must be rotated to 
describe the potential of Eq. |^ When e > 1 and Uq > 
0, the potential is trapping in the zy plane and anti- 
trapping in x; and for e < and Uq > 0, the potential is 
trapping in z and x, but anti-trapping in y. Whereas, at 
e = and e = 1, the potential is independent of x and y 
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respectively. The frequency along the z-axis is 



2eUo 



(21) 



The frequencies along the x and y axis are discussed 
in connection with the net linear microtrap potential 
below. 



^.2 Linear Microtrap Static Potential Analysis 

The static potential is computed using a three-dimensional 
finite element solver. The distance from the center of the 
trap to the bounding box that was used in the simulation 
was more then twice the tip-to-tip cantilever separation. 
To reduce the error in the simulation results, several dif- 
ferent grids were used and the results were averaged. 

It is possible to approximate the three-dimensional 
static potential of the linear microtrap, U-^-^,^{x, y, z) by 
doing a Taylor expansion about the center of the trap. 
Because equal voltages are applied to all capping elec- 
trodes as shown in Fig.[21[b), the cross-terms in the Tay- 
lor expansion are zero. The coefficients of the harmonic 
terms are: 



Dy = 



D, = 



1 




C/o 


dx^ 


1 








1 




C/o 





(0,0,0) 
(0,0,0) 
(0,0,0) 



(22) 
(23) 
(24) 



The derivatives are then evaluated numerically on the 
calculated potential along the axes. The potential is there- 
fore 



f/, 



LMT ^ 2 



D^x'' + Dyy'' + D,z'') 



(25) 



A static potential geometric efficiency factor k compares 
the static potential of the linear microtrap with the hy- 
perbolic electrode geometry of similar characteristic di- 
mension. The characteristic dimension of the linear mi- 
crotrap that corresponds to the distance s in the hyper- 
bolic electrode geometry is the distance from the center 
of the trap to the nearest point on the end-cap elec- 



trodes: des 



KV2 + g)2. 



(26) 



The static potential in the linear microtrap can then 
be written in the same form as the potential in the hy- 
perbolic electrode geometry. 



Uqk I 

"eft 



(27) 



where e = —Dx/D^ = l+Dy/Dz- Given this approxima- 
tion of the electrostatic potential in the linear microtrap, 
the form of the trap frequency along the z-axis is sim- 
ilar to that of the hyperbolic electrodes (Eq. I21|) with 
the difference being only the static potential geometric 
efficiency factor k 



l2K,eUo 



(28) 



The results characterizing the linear microtrap for k 
and e from the numerical simulations are shown in Fig. 

[El 



5 Net Potential 

The combined static and ponderomotive potentials that 
determine the motion of a ion in the linear microtrap 
are written as a three-dimensional uncoupled harmonic 
oscillator potential: 



+ -j2- {-^^ 



(l-e)y2 + z2) (29) 



where the residual axial ponderomotive potential has 
been neglected. Considering this full potential, the effec- 
tive trapping frequencies consist of the quadrature sum 
of the ponderomotive and the static frequencies. 



'^P.LMT - LMT 



' e2y(2r/2 2eKeUo 



P,LMT 



(30) 



(31) 



e2VJf 7;2 2(1 - e)KeUo 



2to2|22£4 
1 en 



1 2KeUo 



(32) 



Table ^ provides a few examples of the calculation of 
the total trap frequencies given a specific geometry. The 
mass of the ion used in calculating the frequencies was 
i"Cd with an RF frequency of i7r/27r = 50MHz. The 
values for rj, e, and k were taken from Figs. Island 1131 
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Fig. 13 The results for the three-dimensional numerical simulations of the static potential in the linear microtrap. Both the 
anisotropy factor e and the static potential geometric efficiency factor k are shown. The ratio of the layer separation over the 
layer width was fixed at one and the gap separation at two {g—2 from Fig. |2jb)). 



a=40/im 
d=10/im 
w=10/im 
b=100^m 


a = 4 

S = 1 
4ft = 21^m 
deft = 59^m 


7? = 0.7 

e = 3 
K = 0.3 


Vo = 40V 

t/o = 20V 
^p.lmt/271- = 20MHz 
'^^.lmt/Stt = 8.7MHz 


Lu^ lmt/27!- = 13MHz 
'^h.lmt/Stt = 23MHz 
'^^.lmt/Stt = 8.7MHz 


a=40/im 
d=2pm 
w=2^m 
b=100Atm 


Q = 20 
S = 1 
4ft ~ 20/^m 
deft — 58/im 


77 = 0.43 
e = 3.5 
K = 0.26 


V^o = 20V 

Uo = IV 
'^plmt/Stt = 6.7MHz 
c^,'lmt/27i- = 1.8MHz 


^^.lmt/Stt = 5.8MHz 
^!/!LMT/27r = 7.3MHz 
^z.lmt/Stt = 1.8MHz 


a=80/im 
d=2^m 
w=2/im 
b=160^m 


Q = 40 
5 = 1 
4ft ~ 40/^m 
deft = 89/im 


77 = 0.38 
e = 3.2 
K = 0.28 


Vo = 35V 

Uo = 0.9V 
Wplmt/Stt = 2.6MHz 
^^lmt/Stt = 1.2MHz 


= 1.5MHz 

^m'lmt/Stt = 3.2MHz 
cj^lmt/St = 1.2MHz 



Table 1 Sample calculations for trap performance. A trap frequency of Ot/^tt = 50MHz was assumed for a ^^^Cd ion. 



6 Microtrap Principal Axes 

Principal axes are the axes along which it is possible 
to describe the motion of an ion in the total potential 
as a three-dimensional uncoupled harmonic oscillator. 
This means that the motion of the ion along each axis 
is independent of the other two spatial coordinates. The 
equations of motion for an uncoupled harmonic oscillator 
are 

X = —oj^x, etc. (33) 

An uncoupled harmonic oscillator corresponds to a po- 
tential with symmetries along the principal axes. Since 



the RF ponderomotive potential fEa. 116(1 is radially sym- 
metric, the principal axes of a linear ion trap are deter- 
mined by the static potential. The principal axes of an 
ion trap are of concern when considering laser cooling an 
ion in the trap. Laser cooling along all three dimensions 
of motion is possible only if the laser wave vector k/a^er 
has a vector component along all three principal axes. 
The symmetry of the microtrap is such that the z-axis 
is a principal axis, therefore, the axes of concern are in 
the xy plane. It is possible to rotate the principal axes by 
applying different static voltages to the electrodes, which 
give rise to an xy cross-term in the static potential. 
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To find the new principal axes, one can rotate the 
coordinate system via Eq. |^ 



X = X cos a + y sm f 
y — —x' sin 9 + y' cos ( 



(34) 



This rotation can be appHed to the potential with an xy 
cross-term of magnitude A 



U' 

LMT 



-ex 



{1 ~ e)y^ + Xxy + z^) (35) 



to find an angle at which the cross term in the rotated co- 
ordinate system [X'x'y') vanishes. This new coordinate 
system, rotated about the z-axis by an angle 0, now de- 
termines the principal axes of the trap. The angle at 
which the cross-term vanishes is found as a function of 
the coefficient of the cross-term A and the geometric fac- 
tor e. 

tan(20) = (36) 

The static axial potential can then be written in an un- 
coupled form, showing explicitly the new principal axes 
x' and y' . 



en 



- (1 - e')y'' + z') (37) 



where e' = ecos(26') + | sin(26') + sin^ 9. 



(a) 



m 



+q 





Fig. 14 (a)Unrotated static 12 point charge potential, 
shown as cross section in the 2 = plane, (b) By changing 
the charge on four of the eight end-cap points, the principal 
axes rotate, (c) The same result can be achieved by apply- 
ing additional negative charge to the center electrodes. The 
aspect ratio of a = 4/3 was used to illustrate the rotation of 
the axes of symmetry of the potential. 



A simple point charge potential model can be used 
to provide a qualitative idea of how the principal axes 
may be rotated. Twelve charges are fixed at the corners 
of three rectangles as shown in Fig. ll4f a). The positions 
of eight charges of value +q are at (±a/2, ±d/2, ±b) and 
an additional four with charge —q at (±a/2, ±d/2, 0). 
A Taylor expansion of the point charge potential where 
a,d can be written as 



U 

^-^ point — 2 
^0 



(1 - e)y' + z') 



(38) 



where Uq = 2g/(47r£o?'o), = ^Jiajiy -t- (d/2)2, and 
e — {2a^ ~ d^)/{a'^ + d^). If two charges are increased 
from q to q' on either end-cap as in Fig. I14r bl , the prin- 
cipal axes are rotated. Alternatively, one could increase 
the negative charge on two of the four point charges in 
the z = plane. This would correspond to applying a 
negative static potential to two of the center electrodes in 
the linear microtrap and is more effective at rotating the 
principal axes. The potential in the point charge model, 
with the addition of these modified charges, becomes 



C/point = '^^^^ {-ex^ - (1 - e)y' + \xy + 7?) 
(47r£o)ro 

where now, the xy cross term has a coefficient 

ad 



A = 6 
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q + q' + d'^ 



(39) 



(40) 



Substituting Eq. ^U] into the condition for the rotation 
angle (Eq. 1361) , and using the explicit form for e in the 
point charge model, the rotation angle can be expressed 
as a function of the applied charges and the trap aspect 
ratio (a = a/d, the tip-to-tip cantilever separation over 
the layer spacing). 



tan 29 = 



1 - q'/q 2a 
1 + q'/q - I 



(41) 



There are several features of this model that give a qual- 
itative understanding of the rotation of the principal 
axes. First, for a given trap aspect ratio a, by increasing 
the ratio of charges, one can rotate the principal axes 
a fixed amount. However, as the aspect ratio increases, 
the amount of rotation that can be given the principal 
axes by changing the charge ratio is decreased, eventu- 
ally approaching zero. 

The principal axes rotation in the xy plane for the 
linear microtrap as a function of the applied voltage on 
two diagonally opposing center electrodes is shown in 
Fig. El The aspect ratio was fixed at a = 20 and 6 = 
1. The other two center electrodes were held at static 
ground with all eight end-cap electrodes at Uo — IV. 
As discussed above, by applying small voltages to the 
appropriate center electrodes, it is possible to rotate the 
principal axes so that laser cooling is effective. 



12 



M. J. Madsen et al. 




•5 5 

Applied Voltaae (V] 

Fig. 15 The principal axis rotation from the xy axis shown 
in Fig. 13 a) as a function of the applied voltage on two diago- 
nally opposing center electrodes. The other center electrodes 
were held at static ground with the end-cap voltages fixed 
at C/o = IV. The trap dimensions are a = 40/im, d = 2^m, 
w=2/im, with all electrodes having a width of lOO/xm. 



7 Conclusion 

A new design for a microfabricated linear ion trap has 
been discussed. Calculations of the RF ponderomotive 
potential have shown a surprising degree of isotropy near 
the center of the trap, even for very high aspect ratios. 
We find that for high transverse electrode aspect ratios, 
the trap strength approaches I/tt times that of a com- 
parable hyperbolic electrode structure. This may be of 
importance in the design of microtraps in applications 
such as Cavity QED and miniature mass spectrome- 
ters where conventional ion trap designs can not be used. 
Geometric scaling factors for the linear microtrap pro- 
vide an easy comparison between these new trap designs 
and conventional ion traps, facilitating implementation 
in future experiments. 

Further investigations will require actual fabrication 
and experimentation with this new type of trap and in- 
clude an investigation of the patch potentials on the sur- 
faces of the doped semiconductors, the limiting electric 
field, and laser scatter from the small aperture. These 
factors are all technical in nature and should not prohibit 
the future implementation of this novel linear microtrap 
design. 



A Appendix: Analytic Solution of the 
Transverse Potential 

The cross-section of the linear microtrap can also be 
modeled as semi-infinite electrodes in a complex plane. 
This model enables calculation of an analytic solution for 
the geometric factor rj in the limit of infinitely thin elec- 
trodes. Following the analysis of parallel-plate capacitor 
fringe fields of Valluri et.al the cross-section of the 
left cantilever electrodes are described in the complex 
plane as lines that go from negative infinity along the 
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Fig. 16 The linear microtrap model in the complex w plane 
with semi-infinite electrodes that terminate at ±a/2 ± zd/2 
with applied voltages ±Vo/2. 



real axis and terminate at — a/2±Z(i/2, as shown in Fig. 
1161 The right set of electrodes (not shown) are a mirror 
image across the x = line and terminate at a/2 ± ic?/2. 
The electrodes are then mapped to an infinite parallel 
plate capacitor. The function that does this mapping is 



± 



2w7r 



ATT 

T 



1 



(42) 



The positive value maps the parallel plate capacitor to 
the left set of electrodes in the w plane, and the nega- 
tive corresponds to the right set. The potential in the 
strip between the two electrodes in the z plane is simply 
the potential between two parallel plates in a capacitor, 
written in complex notation: 



2tt 



Im(z), 



(43) 



where Im(z) denotes the imaginary part of z. 

To find the potential of the original electrode geom- 
etry, the inverse function of Eq.021is needed. With that 
inverse map the potential in the w plane can be eval- 
uated. The inverse map can be written in terms of the 
Lambert W function, Wfc(^), following |25|. 



^± = C± - Wfe (e' 



(44) 



where C± = i^x^ + IT ^ -'^ ^ scaled complex variable. 
The Lambert W function y — Wk{x) is the solution 
to the equation x = yexpy. For complex variables, it 
is important to select the proper branch of Wfe(^) when 
evaluating the function. The appropriate branch is found 
using pS) 

"Im(C) - tt" 



k = 



27r 



(45) 
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where [] denotes the ceihng function which indicates 
that the argument inside the ceihng function should be 
rounded up to the nearest integer. 

If the tip-to-tip cantilever separation a is much greater 
then the layer separation d {a = a/d ^ 1), the poten- 
tial at the center of the trap can be approximated as the 
linear combination of the potential from both the left 
electrodes and the right electrodes 



<..^(lm(.. 



Im(z_)). 



(46) 



With this approximation (a 3> 1) an asymptotic form of 
the Lambert W function exists that leads to a simplifi- 
cation of the inverse map Eq. 1441 The principal branch 
of the Lambert W function has an asymptotic form: 



Wo(0«lne-ln(lnO,e»l- 



(47) 



Inserting Eq. 1471 into Eq. 1441 the inverse map becomes: 
z w lnC±. Expanding the log function about w = 0, Eq. 
^]can be written: 



z± 



±- 



2'KW 



air 



1 


27111] 


2 


2 


air — d 


+ . . . 






( 



(48) 

Since the potential is the linear combination of Im(z+) 
and Im(z_) fEa. l46|l and the linear terms are opposite in 
sign, only the quadratic term contributes to the potential 
of the microtrap. Squaring the complex variable w — u + 
iv and keeping only the second-order imaginary terms, 
one finds that the potential is 



AttVo 
{an — dy 



tUV. 



(49) 



By rotating the coordinate system about the origin by 
9 = 7r/4, the potential is written in a form that allow for 
easy comparison with the quadrupole potential of Eq. 

im 

(50) 



{air — d) 

The geometric factor rj can be found for a microtrap with 
effective distance £cb — \J (a/2)^ -|- {d/2Y. 



47r 



(flTT — d) 
TT 2 

{aiT — 1) 



2 ^cff 



(51) 



This analytic solution of the geometric factor is valid in 
the limit where the trap aspect ratio is large: when the 
tip-to-tip cantilever separation is much larger then the 
layer separation. The geometric factor asymptotically 
approaches 77 = I/tt in this limit. The analytic solution 
fEa. l51|l is shown as the solid line in Fig.|5l Note that this 
complex model assumes infinitely thin electrodes which 
correspond to a large value for the ratio of the layer 
separation to the layer thickness 5 ~ d/w — > 00. The 



values for 77 found via numerical simulations approach 
the analytic solution as d increases and also approach 
the asymptotic value of 77 = l/vr for large a. 

In addition, the analytic model can be used to calcu- 
late the asymptotic values for the ponderomotive poten- 
tial depth and the maximum trap size along the weak 
axis rmax ■ Inserting the asymptotic form of the Lambert 
W function (Eq. \^ directly into the potential (Eq. \^ 
and evaluating the imaginary part, the potential can be 
written directly as a function of u and v. 



2tt 



tan 



2v 



-tan 



2u + a - 

-2v 



d/n 



-2u 



d/n 



(52) 



The pseudopotential can then be directly evaluated, us- 
ing a two-dimensional gradient, from Eq. |H1 The maxi- 
mum of the pseudopotential along the u-axis {v = rmax) 
lies at 



= ^{an-d) 
a , 1 , 



(53) 



The location of the potential maximum asymptotically 
approaches r^ax — a/2 as the aspect ratio goes to infin- 
ity. The trap depth is the pseudopotential evaluated at 
this maximum: 



V'(''n 



1 



4ml22 ^2^2 (1 _ 

V CKTT > 



(54) 



The analytic solution for the scaled trap depth is shown 
in Fig. [71 and approaches the asymptotic value of 2694 
[K-^m^/V^] for large a. 
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